Ginzburg-Landau theory of superconducting surfaces under electric fields 
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A boundary condition for the Ginzburg-Landau wave function at surfaces biased by a strong electric 
field is derived within the de Gennes approach. This condition provides a simple theory of the field 
effect on the critical temperature of superconducting layers. 



The critical temperature of a thin superconducting 
layer is increased or lowered by an electric field applied 
perpendicular to the layer. Similarly to the conductiv- 
ity of inverse layers in semiconductors, superconductivity 
of thin metallic layers can thus be controlled by a gate 
voltage, which makes these structures attractive for ap- 
plications. 

In this paper we show that the phase transition in 
a thin metallic layer is conveniently described by the 
Ginzburg-Landau (GL) theory, where the electric field 
E enters the GL boundary condition as 
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Briefly, the logarithmic derivative of the GL function ip 
or the gap function A at the surface is a sum of the zero- 
field part l/6o and the field induced correction E/Ug- 

The zero-field part has been derived by de Gennes^ 
from the BCS theory. A typical value 6o ~ 1 cm is 
large on the scale of the GL coherence length, therefore 
this contribution is usually neglected. This approxima- 
tion, l/6o « 0, corresponds to the original GL condition 
VV' = 0. 

Here we employ the de Gennes approach to derive the 
field induced correction E/Us- The correction becomes 
important for the above mentioned experiments, where 
fields of the order of 10^ V/cm are realized. Small elec- 
tric fields appearing e.g. in Josephson junctions do not 
require such corrections. 

We start from the condition 
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derived by de Gennes (Eq. (7-62) in Ref. 6). Here No 
is the density of states of a bulk material, V is the BCS 
interaction, and N{x) is the local density of states at 
position X. The actual gap function A{x) has a non- 
trivial profile close to the surface at a; = 0, but it has only 
slow variation at distances exceeding the BCS coherence 
length ^0 — 0.18 hv^ /k^Tc. For x ~ it is crudely linear 
A(a:;) « Aq (1 + x/h), so that Ao is not the true surface 
value but the extrapolation of the gap function towards 



the surface. In Eq. (2) we have used the GL coherence 
length at zero temperature ^(0) = 0.74 for pure metals. 

In measurements of the field effect on the transition 
temperature, the zero-field term 6o is included in the ref- 
erence zero-bias transition temperature. Accordingly, we 
can assume a model of the crystal for which l/6o — 0. 
The simplest model of this kind is a semi-infinite jel- 
lium, where for zero field the density of states is step- 
like, N{x) = No for a; > and N{x) = elsewhere. 
Using that the gap function is restricted to the crystal, 
A(a;) = for a: < 0, one can check that from (2) follows 
l/6o = 0. 

Now we include the electric field. According to the 
Anderson theorem^, the electric field does not change 
the thermodynamical properties directly but only via the 
density of states. The change of the density of states is 
also indirect. The penetrating electric field induces a 
deviation 5n of the electron density. The density devi- 
ation changes the Fermi momentum. Since the density 
of states depends on the Fermi momentum, its value be- 
comes modified. We express this complicated indirect 
effect approximatively via a local linear expansion 



N{x) = No + ^5n{x). 
on 

The de Gennes condition (2) then reads 
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The actual space profile of 5n in superconductors 
is unknown. In fact, some of recent measurements 
suggests that the electric field penetrates deep into 
superconductors.^ Interpretation of these observations is 
not yet settled, therefore we prefer to assume that the 
screening in superconductors is similar to the screen- 
ing in normal metals so that 5n is non-zero only on the 
scale of the Thomas-Fermi screening length. The typical 
Thomas-Fermi length is less then one Angstrom, while 
the gap function varies on a scale typical to the BCS ker- 
nel ~ ^0- Accordingly, in the integral (4) we can take 
A(a;) « A(0) and obtain 
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In this rearrangement we have used the surface charge 
determined by the apphed field eoE = — e dx 5n{x). 

The effective potential Us given by (5) depends on bulk 
material parameters ^o, -^o^ £ind dNo/dn, and on the 
ratio of the gap at the surface to the bulk value 
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According to de Gennes estimates^, the surface ratio is 
of the order of unity. A heuristic derivation of the field- 
effect from the GL equation cannot cover this factor.^ 

It is advantageous to express the effective potential 
Us via the usual parameters of the GL theory. First 
we employ the BCS relation for the critical temperature 
/cbT'c = 1-14 huju exp (— 1/iVoV^). The critical tempera- 
ture depends on the density of electrons. Comparing al- 
loys with different impurity doping, it has been deduced 
that the dominant density dependence enters the criti- 
cal temperature via the density of states. We can thus 
assume duju/dn « and dV/dn « with the help of 
which we express the derivative of the density of states 
via the logarithmic derivative of the critical temperature. 
Formula (5) then simplifies to 
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Here we have expressed the electron density via the Lon- 
don penetration depth A^(0) = rn/(yUone^). Its ratio 
to the GL coherence length defines the GL parameter 

Let us estimate the effective potential Us for niobium. 
The charge carriers are electrons, e = — |e|, with the mass 
close to the electron rest mass, m m 1.2 mc. The GL 
parameter is on the edge of the type-I and II materials, 
K = 0.78, and the logarithmic derivative is of moderate 
amplitude, dlnT^/dlnn = 0.75 (see Ref. 11). Taking 
1] K, \ one finds, Us — —1.3 10^ V. As one can see, a 
large field E ^ 10^ V/cm is necessary to create a field- 
induced correction at least comparable to the commonly 
neglected zero-field value l/&o ~ 1/cm. 

The effective potential (7) is the major result of this 
paper. Now we use it in the boundary condition (1) to 
evaluate the transition temperature T* of a biased layer 
of a finite thickness L. General steps of our analysis 
parallel the theory of the Little-Parks effect^^. It is also in 
a close analogy to the theory of surface superconductivity 
in short coherence length materials^^. 

Let us assume that the electric field is applied only to 
the left surface at cc = 0, while the right surface a,t x = L 
is free of the field. We take l/&o — for simplicity, so 
that we use the boundary conditions 
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The GL function is given by the dimension-less GL 
equation C^(r)VV+V'- 1^1 V = 0, see Ref 12. The tran- 
sition point is characterized by an infinitesimally small 
GL function, jV'P 0. At the transition tempera- 
ture T*, the non-linear term in the GL equation thus 
vanishes ^^(T*)V^?A + V' = 0. This equation is solved 
by ip{x) (X cos [(x — L)/^(T*)], which satisfies the right 
boundary condition (9) while the left boundary condition 
(8) demands 
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When the superconductor has a coherence length ^ 
which satisfies the condition (10), the non-zero GL wave 
function nucleates and the system undergoes a transition 
to the superconducting state. Since the coherence length 
is a function of temperature. 
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one can find from (10) and (11) the transition tempera- 
ture T* . It reads 
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where the function g[T) is a root of ^fg tan ^/g = r. The 
function g is plotted in figure 1. 
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FIG. 1. The dimensionless shift of the transition tempera- 
ture due to the electric field as given by Eq. (12). The exact 
solution of y/ptauy^ = r (thick full line), the linear approx- 
imation g ~ T for thin layers (tangential dashed line), the 
constant approximation g ^ n /A for large suppressive fields 
(dot-dashed line), and the parabolic approximation g ~ — 
for large supportive fields (thin full line). The insert shows a 
detail of the parabohc approximation. 

Although equation (12) is simple by itself, we find it 
useful to discuss its asymptotic solutions. Let us start 
with the experimentally most important limit. The ef- 
fect of the electric field on the transition temperature 
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is rather small and it is most conveniently observed on 
very thin layers. In this case \EL\ <^ \Us \ and one can use 
the linear approximation ^(t) « r shown as the dashed 
line in Fig 1 . Within linear approximation the transition 
temperature (11) simplifies to 
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Formula (13) shows that surface ratio (6) reduces field 
effect on the transition temperature T* . One can com- 
pare (13) with a simple estimate that assumes that the 
induced charge is homogeneously distributed across the 
layer, e5n{x) « —eoE/L. A bulk critical temperature 
modified by the excess charge is then interpreted as the 
transition temperature, T* = Tc + dTc/dn x 5n. Appar- 
ently, the simple estimate does not include the surface 
ratio T]. 

Based on the simple estimate one is tempted to say 
that dTc/dn can be uniquely determined from experi- 
mental data on T*. The present theory shows, however, 
that surface ratio obscures the observed value. 

Among high-T(, superconductors there are many ma- 
terials of large k. From equation (7) one can see that 
these materials have much lower effective potential Us, 
therefore they reveal a much stronger field effect on the 
transition temperature. With these materials it is pos- 
sible to achieve the opposite limit - the regime of thick 
layers \EL\ > |C/,|. 

A measurement in the regime of thick layers has been 
already performed by Matijasevic et al}"^ We will use 
parameters of their sample to evaluate effects expected 
from formula (12). The sample Smo.yCao.aCuaOj, is over- 
doped with Tc reduced to 50 K. The carriers are holes, 
e = |e|, of the mass m = 3.46 me and the density 
n = 5.75 10^'^/cm'^. Based on the authors claim that 
Tc in a monolayer would increase by 10 K at the imposed 
voltage, one can deduce dlnTc/dlnn = —3.12. With 
a typical GL parameter k = 100 one obtains the GL 
coherence length ^(0) = 1.3 nm and an effective poten- 
tial Us = —56.6 V. Note that this is by four orders of 
magnitude smaller than the niobium value. The applied 
field is enhanced by a large dielectric function to an ef- 
fective value E = 7.8 10^ V/cm. For a sample width of 
L = 50 nm and ~ 1 we find \EL/Us\ = 6.9, which 
confirms that these measurements are in the thick layer 
limit. 

In the thick layer limit one has to distinguish whether 
the electric field supports or depresses the transition tem- 
perature. Let us first discuss the depression which ap- 
pears for dTc/dn x e^E/e > 0. In this case EL/Us > 
and the function g approaches the constant asymptotic 
value g 7r^/4 = 2.47 shown as the dash-dotted line in 
Fig. 1. Since g < relation (12) yields the upper 

estimate of the depression of the transition temperature 
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One can sec that the maximal depression is limited by 
the layer thickness. Within the adopted approximations 
the actual value of the electrostatic field does not mat- 
ter once the asymptotic regime is reached. For param- 
eters of Ref. 14 one finds from formula (14) the lower 
estimate Tc - T* < 0.08 K. The formula (12) gives a 
slightly smaller value, as one can guess from Fig. 1. For 
\EL/Us\ = 6.9 the dimensionless shift is ^(6.9) = 1.89, 
i.e., T* - Tc = -0.064 K. Matijasevic et aP^^ reported 
no suppression of the superconductivity. Our estimate 
shows that the suppression is below the sensitivity of 
their method. 

A different situation is met if the direction of the elec- 
tric field is reversed. The electric field then supports 
the superconductivity since dTc/dn x e^E/e < 0. In 
this case EL /Us < and the function g approaches its 
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-T shown as the thin full 



line in Fig. 1. Since g = T^/^^(T*), the coherence length 
^(T*) is imaginary giving the GL function exponentially 
decaying from the biased surface. In this limit, 
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the critical temperature does not depend on the width of 
sample and increases quadratically with the electric field. 
For parameters of Ref. 14 we obtain T* — Tc = 1.6 K in 
a reasonable agreement with the reported shift by 1 K. 

The increased critical temperature (15) is independent 
of the layer width L. This shows that in this limit the su- 
perconductivity above Tc is stimulated in the same way as 
the superconductivity on the surface of an infinite sam- 
ple predicted by Shapiro"'^^. Formula (15) differs from 
Shapiro's formula (9) only by a factor due to the impu- 
rity limited coherence length assumed in Ref. 15. 

We note that the present discussion does not account 
for the charge reservoirs typical to layered CuO materials. 
For a microscopic study devoted to these materials see 
Ref. 16. We also do not assume an eventual effect of the 
electric field on the chemical composition, e.g., due to 
oxygen motion as proposed in Refs. 17, 18. 

In summary, using the de Gennes approach we have de- 
rived the GL boundary condition for a superconducting 
surface exposed to the electric field. This boundary con- 
dition allows one to conveniently evaluate the field effect 
on surface sensitive phenomena from the GL theory. Its 
implementation is demonstrated for the field effect on the 
superconducting phase transition in metallic layers. Our 
approach recovers known features, in particular, that for 
thin layers the transition temperature can be linearly en- 
hanced or suppressed depending on the orientation of the 
applied field. We have found, however, that compared to 
former theories the linear coefficient is modified by the 
value of the gap at the surface. In the limit of thick lay- 
ers we obtain a field induced surface superconductivity 
with the shift of the critical temperature depending on 
the square of the electric field. We also obtain the upper 
limit on the suppression of the critical temperature being 
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independent of the field and inversely proportional to the 
square of the layer width. These features agree with the 
experimental data. 
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